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*vj ■ We investigate the generation of the baryon asymmetry of the Universe within a SUSY SU{5) x T' 

model of flavour, which gives rise to realistic masses and mixing patterns for quarks and leptons. The 
2 I model employs the see-saw mechanism for generation of the light neutrino masses and the baryon 

asymmetry is produced via leptogenesis. We perform detailed calculations of both the CP violating 
lepton asymmetries, originating from the decays of the heavy Majorana neutrinos operative in the 
see-saw mechanism, and of the efficiency factors which account for the lepton asymmetry wash-out 
processes in the Early Universe. The latter are calculated by solving numerically the system of 
Boltzmann equations describing the generation and the evolution of the lepton asymmetries. The 
baryon asymmetry in the model considered is proportional to the Jcp factor, which determines the 

,S^ , magnitude of CP violation effects in the oscillations of flavour neutrinos. The leptogenesis scale can 

^Lc be sufficiently low, allowing to avoid the potential gravitino problem. 
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I. INTRODUCTION 

T-H ' 

JZ^ In the present article we consider the generation of the baryon asymmetry of the Universe in the SUSY model of 

/y-s I flavour based on the SU{5) x T' symmetry, which was developed in [l|, y|- The model possesses a number of appealing 
^^ ■ features which makes it worthwhile to investigate whether it provides also a viable scenario for the baryon asymmetry 
^^ . generation. 

-y-! ' The group T' is the double covering of the symmetry group of the tetrahedron A4 (see, e.g., Q). It was realised 

f~^ , by a number of authors (see, e.g., [j|) that the T' symmetry can be used for the description of masses and mixing of 
p\| ■ both leptons and quarks. The SU{5) x T' model of flavour of interest accounts successfully for the pattern of quark 
masses and mixing, including the CP violation in the quark sector [^. It is free of discrete gauge anomalies [5| and 
gives rise to realistic masses and mixingof the leptons as well. 

The SU{5) x T' model proposed in [l|, y| we will discuss in the present article, includes three right-handed (RH) 
^\ . neutrino fields Nm, I = e,iJ.,T, which possess a Majorana mass term. The light neutrino masses are generated by 
the type I see-saw mechanism and are naturally small. The light neutrino mass spectrum is predicted [6| to be with 
normal ordering and is hierarchical (throughout this article we use the definitions and the conventions given in |7|). 
The neutrino masses rrij, j = 1,2,3, are functions of two real parameters of the model [2|, Q. The latter can be 
determined by using the values of the two neutrino mass squared differences, Am^i and Amg^, or of Am|^ and the 
ratio r = Am2i/Am^i, obtained in the global analyses of the neutrino oscillation data. Using the best fit values of 
Am^i = 7.58 X 10~^ eV^ and Amli ~ 2.35 x 10~^ eV^, found in the analysis performed in [1|, we have p: 

mi = 1.14x lO^^eV, ma = 8.78 x 10"^ eV, mg = 4.867 x 10"^ eV. (1) 

The values of mi, TO2 and m^ are essentially fixed: the uncertainties corresponding to the 3cr ranges of allowed values 
of Am2i and r are remarkably small [6|. 
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The part of the Pontecorvo, Maki, Nakgawa and Sakata (PMNS) neutrino mixing matrix (see Q), resulting from 
the diagonahsation of the Majorana mass term of the left-handed flavour neutrino fields vij^^x), I = e,fj,,T, which 
is generated by the see-saw mechanism, is of the tri-bimaximal form Q- The latter is "corrected" by the unitary 
matrix originating from the diagonahsation of the charged lepton mass matrix Me (for a general discussion of such 
corrections see, e.g., |l(]t4l2| ). Since the model is based on the SU{5) GUT symmetry, the charged lepton mass matrix 
is related to the down-quark mass matrix Md- The model exploits the Georgi-Jarlskog approach for obtaining viable 
relations between the masses of the muon and the s-quark [l|, \^. The Cabibbo angle is given by the "standard" 
expression: 6c = y^md/iris, md and nis being the masses of the d- and s- quarks. As a consequence, in particular, of 
the connection between Mg and Md, the smallest angle in the neutrino mixing matrix ^13, is related to the Cabibbo 
angle (1|: 

sin 6*13 = — = sin6'c. (2) 

Here we implicitly assumed the "standard" parametrisation of the PMNS matrix Q : 

/ C12C13 S12C13 sise""^ \ 

-S12 C23 - C12 S23 siae"^ C12 C23 - S12 S23 sise*'' S23 C13 diag(l, e'^, e'^) , (3) 

\ S12 S23 - C12 C23 si3e"^ -C12 S23 - S12 C23 sise"-^ C23 ci3 / 
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where c^ = cos^i^and Sy — sindij, < 9ij < tt/2, S and aji, j = 2,3, are the Dirac and the two Majorana CP 
violation phases [1J|, < 5 < 2n and, in general |14|, < aji/2 < 2n. We will use this parametrisation in the 
discussion which follows. 

The CP violation, predicted by the model, can entirely be geometrical in origin [2|. This interesting aspect of 
the SU{5) X T' model we will consider is a consequence of one of the special properties of the group T', namely, 
that its group theoretical Clebsch-Gordon (CG) coefhcients are intrinsically complex ^^. The only dominant source 
of CP violation in the lepton sector of the model is the Dirac phase S. The two Majorana phases present in the 
PMNS neutrino mixing matrix, a2i and a^i, are predicted to leading order to have CP conserving values. In the 
standard parametrisation of t/pMNS we have: 021 = and a3i = n. Higher order corrections induce small CP violating 
deviations of the order of few degrees from these CP conserving values of the two phases Q . 

The Dirac phase 6 in the PMNS matrix is induced effectively by the complex CG coefficients of the group T'. 
As we shall see in Section 2, S can take two values in the model considered. One was identified in Q and is equal 
approximately to (5 = 57r/4 = 225°, the precise value being 

5 = 226.9° . (4) 

The second possible value of 6 is given to leading order, as will be discussed in Section 2, by 

6^ -=45°. (5) 

4 

The tri-bimaximal mixing value of the solar neutrino mixing angle 612, which corresponds to sin^ ^12 = 1/3, is 
corrected by a quantity which, as it follows from the general form of such corrections [I0l4l2| ). is determined by the 
angle 613 and the Dirac phase i5: 

1 9 /9 

sin^ 6*12 = - H — sin 6*13 cosS (6) 

o o 

In the SU{5) x T' model considered, ^13 is related to the Cabibbo angle, eq. ([2]). 

The rephasing invariant associated with the Dirac phase 6 [l6| , Jqp , which determines the magnitude of CP violation 
effects in neutrino oscillations [13|, predicted by the model to leading order reads P. UU. Illl| : 

Jcp — — 1= sin 013 sin S = — sin 9c sin S . (7) 

3v 2 18 

For S = 57r/4, eq. ^, the correction to the TBM value of sin^ 9i2 given in eq. ©, is negative and sin^ 9i2 = 0.299, 
where we have used eq. ^ and sin 0c = 0.22. This value lies within the la allowed range, found in the global data 
analysis @. We also have, including the higher order corrections [2|: Jcp — —9.66 x 10~^. If (5 = f , eq. ^, the 
correction to the TBM value of sin^ 9i2 is positive and sin^ 9i2 = 0.37. According to the analyses performed in Q 



and in [1^. the current neutrino oscillation data imply respectively sin^ 6*12 ^ 0.36 and sin^ 6*12 ^ 0.374 at 3a. Thus, 
the case oi S = j is disfavored by the data. For the Jcp factor in this case we get: Jcp = +9.95 x 10^"^. 

Since the neutrino masses, the neutrino mixing angle and the CP violating phases in the PMNS matrix have 
essentially fixed values, the model provides also specific predictions @ for the sum of the three neutrino masses, 

mi + m2 + mg ^ 5.9 x f 0^^ eV , (8) 

as well as for the effective Majorana mass in neutrinoless double-beta decay (see, e.g., (l9l|): 

|<m>| =3.4x 10-3 eV. (9) 

It should be clear from the preceding discussion that the SU{5) x T' model of flavour of interest is remarkably 
predictive: the values of the neutrino masses, the type of the neutrino mass spectrum, the values of the neutrino 
mixing angles and the CP violating phases in the neutrino mixing matrix, as well as the effective Majorana mass in 
neutrinoless double beta decay, obtained in the model are essentially free of ambiguities. The predictions for sin^ig, 
sin 9i2, S and Jcp can be tested directly in the upcoming neutrino oscillation experiments. The value of sin^ig one 
gets in the model, for instance, is relatively small, sin6'i3 = 0.058 ^. It lies outside the 2a, but within the 3a, ranges of 
allowed values of sin6'i3, determined in the global analyses of the current neutrino oscillation data j^, |20|. The results 
of the three reactor De experiments on 613, Double Chooz [21|, RENO [23 and Daya Bay [2^, which are currently 
taking data, can provide a critical test of the model. 

In the present article we investigate the prediction of the SU{5) x T' model of flavour proposed in 111 2| for the 
baryon asymmetry of the Universe. The latter is generated in the model via the leptogenesis mechanism |24l.[25|. The 
dominant source of CP violation in the lepton sector and in leptogenesis is the Dirac phase 5 ^ . Therefore there is a 
direct connection between the baryon asymmetry of the Universe and the CP violation in neutrino oscillations. 

The generation of the baryon asymmetry in the SU{5) x T' model of interest was studied in f27|. However, the 
authors of [27J limited the discussion of the baryon asymmetry generation to the calculation of the CP asymmetries 
in the additive lepton charges, ef, generated in the heavy Majorana neutrino decays, £ = e, ij,,t, i = 1,2,3. They 
based their conclusions on the results obtained for these asymmetries. In the present article we perform a complete 
calculation of the baryon asymmetry, i.e., we calculate not only the asymmetries ej, but also the corresponding 
efficiency factors which account for the effects of the CP asymmetry wash-out processes, taking place in the Early 
Universe. The efficiency factors are computed by solving numerically the Boltzmann equations, which describe the 
evolution of the CP violating asymmetries in the Early Universe. The results we obtain for the lepton asymmetries 
ef do not agree with those found in |27| and our results for the baryon asymmetry contradict the claims made in |27| . 

II. INGREDIENTS 

In the L — R convention in which the neutrino mass terms are written with the RH neutrino fields on the right, the 
superpotential of the model leads P, to the following neutrino Dirac mass matrix, 

2^0 + Vo -Co -Co 

Md = I -Co 2^0 -Co + Vo I CoCo^V = %vu , (10) 

-Co -Co + Vo 2^0 

and to the following Majorana mass matrix of the RH neutrinos, 

1 
Mrr = ( 1 I soA . (11) 

1 

In eqs. ([TU)) and PTj) . £^q, 770, C,o, Co ^iid sq are dimensionless real parameters, A is the scale above which the T' 
symmetry is exact, Y^ is the matrix of the neutrino Yukawa couplings in the basis in which the charged lepton and 
the RH neutrino mass matrices are not diagonal, and v^ is the vacuum expectation value of the "up" Higgs doublet 



^ A larger value of 613 can, in principle, be obtained along the lines discussed in ref. Il2ll . 

^ The Casas-Ibarra matrix l26il . which can be an additional source of CP violation in leptogenesis, is real in the model under discussion 



field of the SUSY extensions of the Standard Model. Thus, the neutrino Dirac mass matrix in the model, Mo^ is real 
and symmetric. As can be easily shown, it is diagonalised by the tri-bimaxinial mixing (TBM) matrix: 

/ 72/3 yvs \ 

Utbm = -v/l76 yi/S -v/172 ■ (12) 

V-v^ V173 7172 / 

We have: 

t^^BM ^'^D Utbm = M'^"^ = diag{3£,o + %, %, 3^o - Vo)CoC'oVu (13) 

where all elements in the diagonal matrix M]iJ°* are real. 

The RH neutrino Majorana mass matrix Man is diagonalised by the unitary matrix S: 

S^ Mrr S^Dn^ diag{Mi,M2, A/3) = sqA diag{l, 1, 1) , Af, > 0, j = 1, 2, 3 , (14) 

where 

/I \ 

5*= I/V2 -i/V2 . (15) 

\0 I/V2 i/V2 J 

and Mj are the masses of the heavy Majorana neutrinos Nj (possessing definite masses), 

N, = S^iNiR + Sji CiNmf = C{N,f , j = 1, 2, 3 , (16) 

C being the charge conjugation matrix. Thus, to leading order, the masses of the three heavy Majorana neutrinos Nj 
coincide, Mj = sqA = M, j = 1, 2, 3. It follows from eq. HH) that S*S^ is a real matrix, so S*S'' = SS"^. 

The effective Majorana mass matrix of the left-handed (LH) flavour neutrinos. Mi,, which is generated by the 
see-saw mechanism, 

M, = -MdM^I^mI , (17) 

{Cl[ = -\i7l A'U Vr + h.c.) is also diagonahzed by the TBM matrix dUl), 

U^bmM.Utbm = diag((3eo + %)', ^o- -(-3Co + %)')^^^^%^ = Q diag{m,,m2, mg) Q^ . (18) 

SqA 

Here Q = i diag(l, 1, ±i) is the matrix which determines, as we shall see, the leading order values of the two Majorana 
phases in the PMNS matrix, and m^ > 0, fc = 1, 2, 3, are the masses of the three light Majorana neutrinos. 



where X = ??o(CoCo''^m) ^-nd Z = Co(CoCo^m)- In what follows we will ignore the overall unphysical factor i in Q. The 
values of rrij given in eq. ((T]) correspond to [0| ^ = ±1.71 x 10~^Vu, and Z = ^7.74 x 10~^u„. 

The charged lepton mass matrix M^ is not diagonal; it is diagonalised by a bi-unitary transformation: M^ = 
VeRM^UJ, where VeR and Ue are unitary matrices and M^ = diag{me,'m^,mT), mi being the mass of the charged 
lepton £, i = e, /i,T. The matrix Ue, which enters into the expression for the PMNS matrix, U = UlUu, diagonalises 
the matrix M^M^. The charged lepton mass matrix A/g (with the corresponding mass term written in the R-L 
convention in terms of the chiral charged lepton fields I'^j^ and Z^^) has the following form [2|: 

/ -(l-*)</)o^^ ^oi^'o\ 

A/e = (l + ^)(/)o^o -3?AoCo 00^0 VdVd^o ■ (20) 

V Co y 

It is related to the down- type quark mass matrix Md via the well-known SU{5) relation: A/g = A/J, with the factor 
(-3) in Me replaced by 1 in Md- The up-type quark mass matrix in the model has the form [^: 

Mu= I (i?W <f>'o' + ll A)^l 2/VoCo I ytvu, , (21) 

yVoCo 1 



In eqs. ([201 and (HH, (po, ipQ, (J>'q, i/jq, yd and yt are real dimensionless parameters, yd and j/t are Yukawa couplings and 
Vd is the vacuum expectation value of the "down" type Higgs doublet of the SUSY extension of the Standard Model. 
Fitting the quark sector observables and charged lepton masses one finds that [2| two of the three angles, present in 
the "standard-like" parametrisation of the matrix Ue, are extremely small, sinSJg = 1.3 x 10^^, sin^js — 1-5 x 10~^, 
while the third satisfies: 

sin 6*12 = X sin 61c. (22) 

o 

Thus, to a very good approximation one can set 9^^ = ^^g = 0, and in this approximation Ue takes the form [6|: 
Ue = $i?i2(6'f2), where $ = diag{l,e"^, 1) and 

cos 6f2 sin 9^2 
Ri2{0t2)= I -sm9f2 cos9f2 | . (23) 

1 

It follows from the above discussion that in the basis in which the charged lepton and the RH neutrino mass matrices 
are diagonal, the matrix of neutrino Yukawa couplings Y^, has the form: 

i; = [/t i; 5 = — [/] Md S . (24) 

Vu 

In the same basis, the Majorana mass term for the LH flavour neutrinos, generated by the see-saw mechanism, is 
given by: 

M, = -vl Y^DJ^^yJ ^UD,U^ , (25) 

where Dj^ = diag{Mi, M2, M3) = Mdiag{l, 1, 1), D^, = diag{mi, m2, m-3) and U is the PMNS matrix, 

U^U^.UtbmQ- (26) 

Using the approximate expression for Ue = <f>i?i2 (6*^2)1 with $ = rfiag(l,e"^, 1) and ^12(^12) given by eq. (P5|) . we 
get: 

/ v/273cf2 + y/T/6sl2e~"f \/l73(cf2 - Si2e"*'^) ^/T/2sl2e-''f \ 
U - v^-^12 - yi76c^2e-'^ VV^ish+che-^n -^fV^c\2e-'^ Q, (27) 

V -yi76 7173 7172 / 

where cf2 = cos 0^2 j ■^12_^ sin 0^2- 

As was shown in [a,ll2|j the phase ip in eq. (|27p and the Dirac phase 5 in eq. Q are related as follows: 

d^ = (^ + TT . (28) 

Comparing the expressions in the left-hand and right-hand sides of the equation M\ Me = Ue[Mf)'^U'l and assuming 
first following [3| that i'^i^'i^'^iiC^i^ < and yVoCo < (with cos 0^2 sin 6*^2 > 0), one finds p|: 

^=J, 5=^n. (29) 

The choice (/'oV'oV'oCo ^ together with the choice yVoCo < (see eq. (|2T|) ) allows to get the best description of the 
quark masses and mixing, possible in the model considered. However, one gets similar description also in the case of 
(/)oV'o'0oCo > and y'ipoCo > '^. In this latter case we get for ip and S: 

p=^±n, <5 = J. (30) 



^ This observation is based on numerical results obtained by M. Spinrath. We thank M. Spinrath for communicating to us the results of 
his numerical analysis. 



Numerically, for kp = 7r/4 and s\^ == 0.22/3 (see eq. ([22])), the PMNS matrix, eq. (H?]), reads: 

/ 0.836e-^i-452° 0.546e'3"9° O.OSlSe-^^s.ooo" \ 
V ^ 0.3676*173.380° o.607e'2-829° „ 0.705 Q. (31) 

V -0.408 0.577 0.707 / 

Taking into account the corrections due to the non-zero values of the angles d\^ and 6*33 in \j\ on finds [2| : 

/ 0.838e-'i-626° 0.543e'3-55i° 0.0582e-*45.ooo" \ 
U ^ 0.362e*"2.463° o.610e'3i60° _ o.705 Q. (32) 

V -0.408 0.577 0.707 / 

Obviously, the differences between the approximate and the "exact" matrices ((3T|) and ([32]) are negligibly small. 

The leading order predictions of the SU{5) x T' model for sin ^13, sin 612 and sin 623 were given in the Introduction 
(see eqs. (HI), ([1]), (O and the related discussions). They can be obtained by comparing eqs. (jH]) and ([?7|) and using 
eq. (Hll). 

Equations pip and p2p allow to determine the values of the Majorana phases a2i and a^i. In the parametrisation 
in which the PMNS matrix is written in eqs. ^7\. ([5T|) and ([5^ they are fixed by the matrix Q = diag{l, l,±i) 
and read a2i/2 = and a3i/2 = tt/2 or 37r/2. Thus, 021 and a3i are CP conserving. Note, however, that the 
parametrisation of the PMNS matrix in eqs. ([?7|) . (^1]) and ([5^ does not coincide with the standard one. Thus, in 
order to get the values of the Dirac and Majorana phases 5 and 021/2 and 0131/2 of the standard parametrisation 
of the PMNS matrix, one has to bring the expressions (pi]) or ([5^ in a form which corresponds to the "standard" 
one in eq. ^. This can be done by using the freedom of multiplying the rows of the PMNS matrix with arbitrary 
phases and by shifting some of the common phases of the columns to a diagonal phase matrix Q. The results for the 
"approximate" and "exact" numerical matrices, eqs. (PT|) and ([5^ . is: 

/ 0.836 0.546 O.OE 

f/~ I -0.367e-*3.48° o.607e'i-38° o.705 \ QaQ, (33) 

and dl 



0.836 


0.546 


0.0518e-'226.69 


-0.367e-*3.48' 


0.607e'i-38° 


0.705 


0.408e*3.i4" 


-0.577e-*i-45° 


0.707 


0.838 


0.543 


0.0582e-'226.93 


-0.362e-*3.99' 


0.610e*i-53° 


0.705 


0_408e*3-55° 


-0.577e-'i-63° 


0.707 


-l)andQe = 


dzaff(e-'3.55° ga.63°^_l)_ ^^^ ^ 



C/~ -0.362e-*3.99 o.610e*i-53 0.705 \ QeQ, (34) 

\ 0.408e*3-55° -0.577e-'i-63° 0.707 / 

where g^ = dm5(e-*3i4°^ e*^'*^", -1) and Qe = diagie-'^-^^^e'^-^^" ,-1). Now comparing eq. ^ and eq. ^ with 
eq. ([3]) we can obtain the "approximate" and "exact" values of the Dirac and the two Majorana phases of the standard 
parametrisation of the PMNS matrix, predicted by the model. For the Dirac phase, for instance, we find, respectively, 
S = 226.7° and [3 <5 — 226.9°. Note that the Majorana phases 021/2 and 031/2 in the standard parametrisation are 
not CP conserving Q: due to the matrix Qa (or Qe) they get small CP violating corrections to the CP conserving 
values and tt/2 or 37r/2. 

The high precision provided by the expression (j27p for the PMNS matrix is more than sufficient for the purposes 
of our investigation and we will use it in our further analysis. This allows to get simple analytic results for the CP 
violating asymmetries, relevant in leptogenesis, which in turn makes transparent and easy to interpret the results we 
are going to obtain. 

Equation ([?5)l . as is well known, allows to express Y^ in terms of U, Di,, D^ and an orthogonal (in general, complex) 
matrix [Hi R, R^ R = RR^ = 1: 

Y^^ —U ^/i\,R^/d^. (35) 

From eqs. (|24p - ([55| and p^ . we obtain the following exact expression for the matrix R: 

R = ( /d:) - 1 Q* M*"^^ [/^s a/ S ( ^/d^) - 1 . (36) 

Using the explicit forms of Q = diag{l, 1, ±i), Mj^""^ , Utbm, S and Dj^ ~ M diag{l, 1, 1) we get: 



R = 



' V3 V3 * 

!_ /¥ 

/3 ys 

V 0-1/ 



^ V3 



(37) 



The same expression for the matrix R was obtained in [23|. Thus, in the SU{5) x T' model considered, the R matrix 
is real, i.e., CP conserving [2^ (see also [20|), and symmetric, R* = R, R^ = R, and the elements Rks = Rak = 0, 
k = 1,2. We note that the signs of the entries in the 1-2 sector of R depend on the signs of X and Z: the signs in eq. 
([57)1 correspond to X > and Z < (see eq. ([TO)) and the related comments). 

III. RADIATIVELY INDUCED LEPTOGENESIS 

As we have seen, the three heavy Majorana neutrinos Nj are degenerate in mass at the scale Mx at which the 
Majorana mass matrix of the RH neutrinos is generated. We will assume that this scale does not exceed the GUT scale, 
Mgut = 2 X 10^^ GeV: Mx < Mqut- Actually, in the SUSY SU{5) x T' model considered, we have Mx = Mqut- 
Given the fact that the R matrix is real and GP conserving, the baryon asymmetry can only be generated in the 
regime of flavoured leptogenesis |3C)l.l3l|. The regimes of 2-flavour and 3-flavour leptogenesis are realised, in general, for 
values of the masses Mj = M, j = 1, 2, 3, of the heavy Majorana neutrinos satisfying J23| M ^T < (l+tan^ /3) x 10^^ 
GeV and M ^ T < (1 + tan^ /3) x 10^ GeV, respectively, where T is the temperature of the Early Universe and 
tan/3 = Vu/vd is the ratio of the vacuum expectation values of the two Higgs doublet fields, present in the SUSY 
theories, v = \Jv\ + v\ ~ YJ'l GeV. If the heavy Majorana neutrinos would be degenerate in mass at the scale 
(temperatures) at which the flavoured leptogenesis can take place, as is well known, no net baryon asymmetry would 
be generated. However, if leptogenesis takes place at a scale Mpj^Q < {■^)Mx, higher order corrections accounted for 
by the renormalisation group (RG) equations describing the change of the masses Mj with the change of the energy 
scale from Mx to Mflg ^ (1 + tan^ /3) x lO'^^ GeV, lift the degeneracy of Nj |33435| . generating relatively small 
splittings between Mi, M2 and M3: AMij(MFLG) = ^/^(A'/flg) - Mj^Mflg) ^ 0, i ^ j = 1,2,3. Since the mass 
splittings |AMij(MFLG)| thus generated are exceedingly small, we expect the baryon asymmetry to be generated in 
the regime of resonant flavoured leptogenesis |28l . |36| . 

In the case of resonant flavoured leptogenesis, the GP violating asymmetry in the lepton charge Li, I = e,iJ.,T, 
generated in the out of equilibrium decays of the heavy Majorana neutrino Nj taking place at the scale MpLc is 
given by [28|: 

, _ r{N, ^ e- H+) + r(iv, ^ i/, ijo) - r(iv, ^ e+H-) - t{n, -^DeH") _ i ^ 



Here 



V{N,^ i- H+)^T{Ni^vtH°)+V{Ni^ l+H-)+T{N,^vtH°) Stt ^ ■''' 



Im[{Y^Y,),,{Y,)* {¥,),. ^ 



-H+) 


+ T{Ni -^ VI H°) + 


T{N,- 


'? - 


M.MjAMf^ 


T^ 


Olj 


{AMf^)^ + MfT] ' 


^i] 



{YJY,) 



(39) 



where Y^ is defined in eqs. 



and 



r. = i-(n^>^)«M., (40) 



Stt^ 



t2 _ ,,r2 ,,r2 ^ r,7,,.r2rAf rW _ ^^0 



AMl,^Mf-Mt9^2MfS^, 5-^^-1, J^^. (41) 

The parameter (5^ describes the deviation from complete degeneracy of the masses of the the heavy Majorana neutrinos 
Nj and Ni. All quantities which appear in eqs. (|38)) - (|4T|) should be evaluated at the leptogenesis scale Mflg- The 
baryon asymmetry is generated in the regime of resonant leptogenesis if at Mp^Q the following condition is fulfilled: 

Af,r, -AM|,, z^j. (42) 

We have discussed above the asymmetry generated in the decays of the heavy Majorana neutrinos A^,; into the Higgs 
and lepton doublets. A lepton flavour asymmetry e\ is also generated from the out-of-equilibrium decays of Ni in the 
Higgsino and slepton doublets I. Similarly, the sneutrinos Ni generate CP asymmetries e^ and e~ with, respectively, 

£ and £ in the final state. As can be shown, one has neglecting soft SUSY breaking terms: el — tl — t~ — e~. 

It follows from eq. (|39p that the necessary conditions for a successful resonant flavoured leptogenesis include: i) 
the presence of CP violating phases in the matrix of neutrino Yukawa couplings Y^; ii) non- vanishing off-diagonal 
elements of the matrix Y^Y^: {Y^Yy)ij ^ for z 7^ j; iii) non-degeneracy of the heavy Majorana neutrino masses M^: 



(5^ 7^ 0, i 7^ j. The first requirement is fulfilled by the presence of the CP violating phases in the neutrino mixing 
matrix U. The second and third general requirements are satisfied, as we are going to discuss next, owing to the RG 
corrections in the quantities Mi and Y,y, which have to be included when the latter are evaluated at the leptogenesis 
scale MpLG- 

The RG running of the heavy Majorana neutrino masses Mi depends on the quantity Y^Yi, [33 ■ It proves convenient 
to work at the scale Mx in a basis of the heavy Majorana neutrino fields in which the matrix Y^Yn is diagonal. This 
can be achieved by performing an orthogonal transformation of Nj. The latter can be done without affecting the 
heavy Majorana neutrino mass term since at the scale of interest the heavy Majorana neutrinos Nj are degenerate in 
mass. The change of basis, Nj — Oji^Nj,, where O is an orthogonal matrix, implies the following change of the matrix 
of neutrino Yukawa couplings: Y^ = Y^O. Using eq. (|35p and the facts that D^ = M diag{l, 1, 1) and the matrix 
R is real and orthogonal, there always exists an orthogonal matrix O such that RO, and correspondingly Yl'^Yl, are 
diagonal matrices. Taking into account the explicit form of the matrix R in the model considered, eq. (j37p . in what 
follows we will use 



O 



cos 9 sin 6 

— sin 6 cos 6 

1 



/ ,/I J- o\ 



V 







(43) 



1/ 



It is easy to verify that RO = diag{—l, 1, —1) = /, Ijj 
couplings yj is given by: 



= T]j, with ri2 = — ??i.3 = 1. The matrix of neutrino Yukawa 



Y' = Y,n = 




(44) 



In the new basis we have {Yl^YDij = for i 7^ 7, and 



(>;'^n). 



M, 



■mi 



r' 



Mf 
8^ 



■m. 



M, = M . 



(45) 



where F^ is the N^ total decay width. 

The expression for the CP violating asymmetry ef in the new basis in which Y^^Y^ is diagonal at Mx can be 
obtained from eqs. ([55|) - PO)) by replacing Y^, and F^ with Y^ and F^, respectively. Note, however, that in the new 



basis we have Xf ■ = 0. Thus, the CP violating asymmetries ef will be zero unless non-diagonal elements of YJ^Y^ are 
radiatively generated at the leptogenesis scale Mpj^Q < {<^)Mx- 

As will be shown later, in the model considered a non-zero baryon asymmetry can be produced only in the regime 
of 3-flavoured leptogenesis, i.e. for M < (1 + tan^ /3) x 10^ GeV ^ 4.9 x 10^^ GeV, where we have used the constraint 
tan^ ^ 70 (see '', e.g., [s^l). Taking into account that m, ^ 5 x 10"^ eV and v = 174 GeV, we get |(rjty-j)..| ^ g x 
10-3 < 1. 

In the new basis, the running of the heavy Majorana neutrino masses Mi is governed by the following equation |33| : 



dM, 
dt 



iiYl^Y^uM, 



t 



?r In 

167r2 Mx 



(46) 



where the initial conditions are at the scale // = Mx at which Mi = M , i = 1,2, 3, and the masses Mi are evaluated at 
the scale jjl = Mflg < ('C)A/x. The latter coincides, up to negligibly small corrections, with M: Mflg == M. The 
running of the masses Mi from Mx to Mflg — M induces the splitting between the masses of the heavy Majorana 
neutrinos, necessary for a potentially successful leptogenesis. The solutions of the equations (|46|) |34l . l35l | lead for 
|(Fjtyj)^^| <^ 1 to the following expression for the mass splitting parameter (5^: 






-^[{Yl^Yl),,-{Yl^Y,%]V 



M 



-{ruj -mi)t, j y^i, 



' iJ^'- 



167r2 V M 



(47) 



For Mx/M == 2 X 10^ 2 x 10^ and 2 x lO'^, we get t = 0.092, 0.077, 0.063. The corresponding values of S^, are given 
in Table U 



In the calculation of the baryon asymmetry we will values of tan/3 ~ 10, which are much smaller than the quoted maximal value. 





Mx/M = 2 X 10'' 


Mx/M = 2 X 10^ 


Mx/M = 2 X 10* 


021 


-9.28 xlO"'^ 


-7.81 xlO"" 


-6.33x10"^ 


O31 


-5.80x10"'^ 


-4.88 xlO"^ 


-3.96 xlO"* 


O32 


-4.87x10"'' 


-4.10 xlO"^ 


-3.33 xlO-" 



Table I: Values of the heavy Majorana mass splitting parameter Sl 



The elements of the matrix of neutrino Yukawa couphngs Y/, also evolve with the scale fi when the latter diminishes 
from Mx to Mflg — M. This change is governed by the RG equations for {YD^i, whose general form was given in 
[H-llll. In the case considered by us we have at Mx- {Y/^YD^j = 0, z 7^ j, and ''' |(yj^yj)M| ^ 8 x 10^^. In this case 
the RG equations for (FJ)fi [35| simplify considerably and read: 



dt 



E Vl-l9l-'^9l 



-vl 



{yi\ 



H , 



e,fi,T, 1 = 1,2,3, 



(48) 



where y^, q = u,c,t, and y^, i = e,fi.,T, are the charge 2/3 quark and charged lepton Yukawa couplings, 91^2 are 
the U{1)y and SU{2) gauge couplings of the Standard Model and we have neglected terms ex Y/Z^Y/,. The quantities 
which appear in the square brackets in the r.h.s. of eq. (pS]) evolve with the scale jj, as it decreases from Mx, but the 
effects of their evolution are subdominant for the problem under study and we will neglect them. Thus, we will use 
their values at the scale Mx, which will be assumed to be close, or equal, to Mqut- 

We are interested in the quantities {Yl'^Y/,)ij, i -^ j, at the scale Mflg — ^, which enter into the expression for 
the CP violating asymmetry e^. Since these quantities are zero at Mx, they can get non-zero values at Mflg due 
only to the term involving the charged lepton Yukawa coupling yj in the RG equation (|^5)) |34l . l35j . The solutions of 
the RG equations (|48p in the leading logarithmic approximation lead to the following result: 



(>;'^K).. = -2y?(n'*)r.K)r,t 



^J- 



(49) 



vi = 174 



where yr — {rrir/vd) = {mr/v)^\ + tan^ /3 is the r Yukawa coupling, nir being the r mass, and v = \/v^ 
GeV. Neglecting relatively small effects, the quantities in the r.h.s. of eq. (^^ can be taken at the scale Mx- Note 
that even though at Mx the off-diagonal elements of Y/^Y/, are zero, they have non-zero values at the leptogenesis 
scale MpLG due to the radiative corrections. 

Using the result obtained for {Y/^'^Yl)ij, eq. ([49l), and eqs. ([38]), (|44]) and ([45l), we get for the CP violating 
asymmetry: 



1 

Stt 



vlt 



E 



{S^Y 



s:N 



IGttv'?, 



Mjrrij 



im [u:. 



Urj Ul 



Uij] 



(50) 



It follows from the expression ([50)) for e^ we have derived that e^ +ef + e[ = 0, i = 1, 2, 3. This result is a consequence 
of the fact that the R matrix in the model considered is CP conserving (see, e.g., [28|). One can easily convince oneself 
using the explicit expression for the PMNS matrix p7)) that we also have: e[ = 0, i = 1,2,3. The same conclusion 



is reached also if one uses the PMNS matrix in which the higher order corrections have been included °, eq. ([32| or 
([34| . Thus, in the SUSY SU{5) x T' model of interest the baryon asymmetry can be generated only in the regime of 
3-flavoured leptogenesis |31| . 



^ For a matrix of neutrino Yukawa couplings Yi, such that Re[{YjY,^)ij] 7^ 0, i 7^ j, the RG equations for Y,^ have a singularity in the 
case of degenerate in mass heavy Majorana neutrinos l32ll34l [35ll . As a consequence, the quantity {YjY,j)ij, i ^ j, that enters into the 
expression for the CP violating asymmetry e|, does not vary continuously with the scale when the latter changes from Mx to Mp^c- 
This fact was not taken into account in the calculation of the asymmetries ef performed in [27|| . Since in the basis in which we work 

1, 2, 3, the indicated problem does not appear when we consider the RG 

we obtain. The latter is not difficult to verify. 



: 0, i 7^ j, at Mx at which Mj = M, j 



we have {Y/,''Y/,),j 

evolution of (yj)« and of (yj^yj).^.. 

It is claimed in [27| that ej ^ 0, which does not correspond to the result ej" 
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The requirement that the baryon asymmetry is generated in the 3-flavoured thermal leptogenesis regime combined 
with the upper hmit on tan/3 imphes: M fC 4.9 x 10^^ GeV. As is not difhcuh to show, we have for Af pQ 10^'^ GeV: 



(S 



N\2 
ji ) 



> 



1 MjlTlj 

16tt V?, 



(51) 



For M = 10^'^ GeV, (<5^)^ is bigger by a factor of 10 than the term in the right-hand side of the above inequahty. 
Neglecting the correction due to the latter, we get a rather simple expression for the asymmetry e|: 



y'r 



327r ^-^ rrii 



E' 



-Im[U;,Ur,U^,U,,] 



(52) 



where we have used eqs. (|47)) and (jSOj) . 

Expression ([5^ for ef does not depend explicitly on the masses of the heavy Majorana neutrinos and on the RG 
factor t. Thus, the CP-asymmetries ef are entirely determined by the r Yukawa coupling and the low-energy neutrino 
mixing parameters, i.e., the neutrino masses, the neutrino mixing angles and CP violating phases in the neutrino 
mixing matrix. They depend weakly on scales Mx and M, e.g., via the running of the t Yukawa coupling. The 
asymmetries ef depend quadratically on the r Yukawa coupling and thus on tan^ /3. This dependence is crucial for 
having a viable thermal leptogenesis in the SU{5) x T' model of flavour under consideration. 

From eq. ((52| . using eqs. p7)) and ([7]), we obtain: 



y2 

327r -^ TO, 



j¥'' 



TO, 



Pj^ 



1,2,3, 



(53) 



where pji = —pij, i ^ j, and p2i = P3i = P23 = +1- Thus, the CP violating asymmetries e^'^, i = 1,2,3, are 
all proportional to the Jcp factor, which determines the magnitude of CP violation effects in the flavour neutrino 
oscillations. 

The final lepton number asymmetry, which is partially converted into a non-zero baryon number asymmetry by 
the fast sphaleron interactions in the thermal bath in the Early Universe, receives a contribution from the out- 
of-equilibrium decays of the three heavy Majorana neutrinos (sneutrinos) A^/ (N'i), which are quasi-degenerate in 
mass. The amount of matter-antimatter asymmetry predicted by the model is computed numerically by solving 
the corresponding system of Boltzmann equations. We report below the relevant set of Boltzmann equations in 
supersymmetric leptogenesis [3a, [3^ for the lepton fiavour (lepton charge) asymmetries Ya^ = Ya^ + ^a~, with 



^e{i) 



B/3 



'i{ty 



dY, 



N' 



dz 

dYwu 
dz 

dYA, 
dz 



Y, 



N' 



sH{Mflg) 



sH{Mflg) 



2 {ih + 7s,AL=l) I ^cq 



N' 



^hh + lk 



AL=1 



Y; 



N'i 



Y~ 

^ N'i 



- 1 



sH{Mflg) ^ 

/ i,l , ii 
16+16 



(4 + 4) (71, + 7^, A..i) ( ^ - 1 ) + (4 + 4) {16 + Ik 



lw,AL=l + 7w AL=1 + 



N[ 

16 + Id 



i,l I i,£ 

lw,AL=l + lw,AL=l 





(54) 




(55) 


L = lj T^cq 
\ NU 


\ 

1 
/ 


Er Aw Ya„ ' 


.(56) 


Vcq 



Here Y^/ (Y^) is the N[ (A^^'-equilibrium) abundance, z = Mfi^q/T, T being the temperature of the thermal bath, s 

is the entropy density and H{T) is the expansion rate of the Universe. The quantity 7)^ {i = 1, 2, 3) is the thermally 
averaged total decay rate of the Majorana neutrino N^ into the SM lepton and Higgs doublets. Similarly, 7^ ^j^^i 
is the corresponding AL = 1 thermal scattering rate of N^ with SM leptons, quarks and gauge bosons. The flavour 



"^ As was pointed out earlier, the CP asymmetries eT (i = 1, 2, 3) are equal to zero in the model we are discussing. Nonetheless, a source 
term for ^T(f) is provided by non-zero Y^^ , as is explicit from the flavoured Boltzmann equation 11561 . 
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M/T 



Figure 1: Solution of tlie Boltzmann equations (|54p - (|56p for tan/3 = 10 and 5 — 7r/4. See the text for details. 



dependent washout processes involving N^ inverse decays and the relative AL = 1 scatterings are denoted as 7^' 

and "fw AL=i^ respectively. Finally, the matrix elements of A in supersymmetric type I see-saw scenarios are [3 
Aap = i6/2133 for a ^ /3 and A^a = -221/2133 (a = e, fi, r). 

The entropy density, s, and the expansion rate of the Universe, H{T), are given by: 



i,t(i) 



9 



,27r2r3 



45 



H{T) 



1.66 ^T 



mpi 



(57) 



where g^ = 228.75 Hil and mpi ~ 1.22 x 10^^ GeV is the Planck mass. 

In the case in which the soft SUSY breaking terms are negligible, the thermal rates in ([54|) - ([56)) satisfy the conditions 

|3l| : Jx = 7V ^11*^ 7!y = 7x ~ 7x — 7x • ^^ ^ good approximation, supersymmetric leptogenesis proceeds as a 
manifest generalization of the standard leptogenesis scenario of the type I see-saw extension of the SM. Indeed, new 
effects due to different supersymmetric equilibration mechanisms between particle and sparticle number densities 
provide typically only relatively small corrections [33 , which can be safely neglected for the purposes of the present 
study. 

The dominant contribution to the production and damping of the lepton asymmetries is generally provided by 
decays and inverse decays of N- [40| , whose thermal averaged rates are 



Id 






Id 



Id 



l(n')£ 



{Yl^YD 



(58) 



where lCi{z) is a modified Bessel function of the second kind. 

We neglect in ((54|) - ([56|) . for simplicity, thermal corrections to the CP asymmetries and the decay/scattering rates 
|40| . We do not include either the Ai — 2 washout of the flavour lepton asymmetries in the Boltzmann equations 
listed above because they are subdominant at the temperatures at which the 3-flavoured leptogenesis takes place. ® 

The final baryon number density (normalized to the entropy density of the Universe) is: 



rs-f,^r.. 



Ya 



Ya. 



(59) 



In order to have successful leptogenesis, the CP asymmetries ef {i = e, n) should be sufficiently large and should have 
the correct sign. According to eq. (|53p . the sign of ef = — e^ and, consequently, of Yg, is fixed by the value of the 



As is well known, AL = 2 scatterings mediated by N( {N',) can be safely neglected if r'-/H{T) <C 10 X Mi/ (iQl* GeV) [H 
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rephasing invariant associated to the Dirac phase 5, Jcp. Numerically, from ([55]) . we get for tan^ (3^1: 

e^ ~ -2.3 X 10"^ Jcp (tan /3)^ e^ ~ 1.3 x 10"^ Jcp (tan ^)2 , e;^ ~ 2.1 x 10"^ Jcp (tan /3)2 , (60) 

where we have used eq. ([T]) and y'^ ~ 10~"'tan^ (3. Taking, more explicitly, tan/3 = 10, one easily obtains: 

e^ ~ -sgn(sin (5)1.4 X 10"^ e^ ~ sgn(sin(5) 7.0 x 10"^ e^ ~ sgn(sin^) 1.3 x 10"^ (61) 

which, in general, is the right order of magnitude of the CP asymmetry in order to have a successful leptogenesis. 
Notice that sgn(sin(5) is equal either to (—1) or to (+1), depending on whether the Dirac phase S ^ 57r/4 or S = 7r/4, 
which are the two approximate values S can have in the model considered (see eqs. ([29)) and (pO| ). 
Taking into account eq. (|60p . expression ([5^ can be recast in the form: 

Yb « Jcp {tsmPferjB Y^"} (z < 1) , (62) 

where Y^ = 45/(7r'*5,) ~ 2 x 10"'^, e = 10^^ and i]b > is, by definition, the efhciency factor of the asymmetry. It 
follows from eqs. ([7]) and (|62|) that for S = 226.93° ~ 57r/4, the baryon asymmetry has the wrong sign. Thus, the 
observed value of the baryon asymmetry can be obtained in the model considered only for S ~ 7r/4. 

The the efficiency factor rjB in eq. ([62| can be computed by solving the full system of Boltzmann equations ([54]) - ([56)) . 
We note that in the model considered the parameter i]b does not depend on the leptogenesis scale Mflg ^ M. This 
can be easily understood if one considers, for simplicity, the solution of the Boltzmann equations where only decay 
and inverse decay processes are included: as we have already mentioned, this is a good approximation in thermal 
flavoured leptogenesis. In this case, from eqs ([45]) . ([57| and ([58| one has: 

Z'yl, m.mpi 

oc — . (63) 



sH{M) 



Therefore, the Boltzmann equations do not explicitly depend on the heavy Majorana neutrino mass scale M within 
the indicated approximation. We verified numerically that the dependance of -qs and Yb on M is relatively weak 
also if we take into account the scattering processes. This implies that, in the class of SUSY see-saw models of the 
type considered in this paper, the leptogenesis scale Mpi^Q can be lowered sufficiently in order to avoid the potential 
gravitino problem ^ . 

In figure[TJ we report the solution of the full set of Boltzmann equations (|54p - ([55|) for tan/? = 10 and 5 = 7r/4. The 
red, blue, green and black lines represent |i^^|, \Ya |, I^a^I and \Yb\, respectively. The dashed line corresponds to 
y^, while the other three black lines are the RH neutrino abundances Y/v' ■ The gray horizontal band gives the 3a 

interval of experimental values of Yb'- Yg^'^ = (8.77 ± 0.21) x 10~^^ [42|, where we have quoted the la error. In this 
numerical example, we get the final asymmetries: 

Ya, ^ 4.7x10"^", Ya^ ^ -5.8x10"", Fa, ^ 2.6x10"" and Yb ~ 1.4x10"^°. (64) 

From eq. (j62p and the numerical value of Yb thus computed, we get an efficiency factor tjb — 0.07. Obviously, one 
can get a value of Yb closer to the mean best fit value Y§^^ ~ 8.77 x 10~^^ for a somewhat smaller value of tan/3. 

We would like to conclude with the following remarks. As we have shown, the correct sign of the baryon asymmetry 
in the SU{5) x T' x Z\2 x -^12 model considered [l|, [3 can be obtained only in the case oi 8 = 7r/4. As has already 
been discussed in the Introduction, for this value of the Dirac phase b we have sin^ 6*12 ^ 0.37, while the current 
neutrino oscillation data imply at 3cr sin^ ^12 ^ 0.36 Q, or sin^ 6*12 ^ 0.374 [I^, depending on the details of the 
analysis. For (5 = 57r/4, the value of sin^ Qyi — 0.299 predicted by the model lies within the \a interval of values 
suggested by the data, but the predicted baryon asymmetry of the Universe has the wrong sign ^° (see eq. (|62p ). 
If sin^ 012 = 0.37 would be definitely excluded by future data, one would have to modify the SU{5) x T' model of 
flavour we have considered in the present article. One possible "minimaf modiflcation could be to lift the degeneracy 
in mass of the the three heavy Majorana neutrinos (sneutrinos) at the scale Mx , at which the flavour symmetry is 
spontaneously broken. This could be achieved, e.g., by replacing the chiral superfield S in the SU{5) x T' x Z12 x Z[2 
invariant superpotential of [2| with a new chiral supermultiplet x, which is a Standard Model singlet and is charged 
only under the discrete group ^(2, with charge w^. The model, therefore, has the same gauge and flavour symmetry 
groups and the same number of flelds as the one discussed in [2|. In this new scenario, the flavour structure of the 



® The Davidson-Ibarra bound |4l| does not apply in the radiative leptogenesis scenario discussed by us. 
^^ Our result for the sign of the baryon asymmetry in the case of (5 = 226.93° ~ 57r/4 contradicts the claim made in [27 
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superpotential naturally generates a Majorana mass matrix (term) for the heavy RH neutrinos at the scale Mx- 
The latter is still diagonalised by the tri-bimaximal mixing matrix Utbau but has non-degenerate eigenvalues. The 
low energy phenomenology, as well as the generation of the baryon asymmetry of this class of models is therefore 
worthwhile investigating, but such an investigation lies outside the scope of the present work. 
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